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Abstract 

In a recent paper a model of quark and lepton masses was pro- 
posed. Without any family symmetries almost all the qualitative and 
quantitative features of the quark and lepton masses and Kobayashi- 
Maskawa mixing angles are explained, primarily as consequences of 
various aspects of 50(10) symmetry. Here the model is discussed in 
much greater detail. The threefold mass hierarchy as well as the rela- 
tions ~ m^, ~ 3m^, m-g ~ "i«/''^f tan^c" — 

\frn^[f^s^ ^cb \J nnPs I "^t ^-^^ Vuh ~ KisV'cfe follow from a simple 
Yukawa structure at the unification scale. The model also gives def- 
inite predictions for tan /3, the neutrino mixing angles, and proton 
decay branching ratios. The (f^ — v^-) mixing angle is typically large, 
tan/3 is close to either m^/m^ or m^/m^, and proton decay is in the 
observable range, but there is a group theoretical suppression factor 
in the rate. 
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1 Introduction 



In a recent paper we proposed a model of quark and lepton masses based 
upon the gauge group SO {10). In that paper it was shown how most of 
the features of the fermion masses and mixings can be understood to be 
consequences of various aspects of SO (10) invariance rather than of some 
family symmetry or arbitrarily imposed texture. 

In this paper we go into greater detail. In Section 2 the basic ideas of 
the model are explained and its structure is reviewed in detail. The results 
of numerical fits are presented in Section 3. Section 4 deals with the predic- 
tions of the model for tan (3 and the neutrino mixing angles. Predictions for 
proton-decay branching ratios are also briefly discussed there (for details see 
Ref. (2)). In Section 5 various technical issues relating to the Higgs sector, 
symmetry breaking, and the naturalness of the gauge hierarchy are exam- 
ined. Finally, in Section 6 certain alternative possibilities for variant models 
are discussed. Various technical details are summarized in three Appendices. 

2 Review of the Model 

(a) The Root Model 

The model we shall study is based on SO (10), and though supersymmetry 
is not essential for its account of quark and lepton masses it shall be assumed 
because of the gauge hierarchy. The "matter" fields (quarks and leptons) are 
contained in three spinors (16j, i = 1,2,3), which are the "families" , and a real 
(in the group theory sense) set of additional representations. These latter 
consist of a pair of family and anti-family (16 -|- 16), and a pair of vectors 
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(10 + 10'). In the "long version" of the model there is also a pair of adjoints 
(45 + 45'). The Higgs multiplets which couple to the matter are in a vector 
(IOh): two adjoints {45h + 45i7), and a pair of spinors {16h + IQh)- The 
short version of the model has two sets of Yukawa terms in the superpotential, 

W^pinor and WvectoT- The long version has an additional set of terms, M4,djoint- 
These are given by 

W^spinor = MT616 + ^6^(16^16)45^ + ^0^(1646)10^^, (1) 

i i 

W^vector = d{10 10')45^/ + ^Q(16aO)16H + X^4(16a0')16^, (2) 

i i 

W^adjoint = /(45 45')45j, + ^e,(16,45)T6H + ^e[{16i45')T6H. (3) 

i i 

Note that all three pieces of the Yukawa superpotential have the same ba- 
sic structure. In each there is one term (the first) which couples together the 
pair of extra representations and gives them superheavy mass, and two terms 
(the second and third) which couple these extra representations to the ordi- 
nary families. What is not allowed is a direct mass term g(y(16il6j)10ij, 
which would, unless there were fine tuning, tend to give comparable masses 
to all the generations. All these three pieces, in other words, have a kind of 
"see-saw" structure in which the ordinary families get masses through their 
mixing with the extra fields. 

The dominant contribution to the fcrmion mass matrices is assumed to 
come from H^spinor- (As will be seen in Section 5, this is most simply explained 
as being a consequence of the condition (5(10ij)) <^ (5(16^^)). Here and 
throughout the notation p(q) refers to a p of SU{5) contained in a q of 
5*0(10).) Diagrammatically, this contribution comes from the tree graph 
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shown in Fig. 1. Its approximate form can be read off from that figure. 
Wo = }_^aibj — (16il6^). 



(4) 



Defining di = ai/a, hi = bi/b (a = \a\,b = \b\); writing the VEV of the 45^^- 
as {4:5h) = ^Q, where Q, hke M, is of order Mqut, and where Q is a hnear 



combination of 5*0 (10) generators; and defining T = one has 



Wo ^ aTilOn) E l^ihQ( 



'16,16o 



(5) 



Then for charge- (|) quarks the mass matrix, Uij is given by 



Uij = aTvY^ [dibjQu + djbiQu- 



(6) 



with similar expressions for the mass matrices of the charge- (^) quarks 
(Dij), charged leptons (Lij), and the Dirac mass matrix of the neutrinos 
(Nij). One can choose, without loss of generality, the axes in family space 
so that hi — (0,0,1) and di — (0, sin^, cos^). The quantity Q is a linear 
combination of 5'O(10) generators, with in general complex coefficients since 
45h is a chiral superfield. There is a two-dimensional space of such generators 
that commute with SU{S)c x SU{2)l x U{1)y- We can thus choose to write 



g = 273R+|e(y/2), 



(7) 



where I^r is the third generator of SU{2)r. Then Eq. (6) and its analogues 
take the form 



Un = aTv 



/ \ 

g„sin^/Ar„ 

\ Q^c sin e/N^o (g„c + g„) COS e/N^cN^ J 



(8) 
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/ \ 

Do^aTv' Qdsm9/Nd , (9) 

V Qdc sin 9/Ndc {Qdc + Qa) cos 9/NdoNd j 



'0 — 



aTv' 




(10) 



A^n = aTv 



\ 

sin 9/N^ 

Q^c sin 9/N^c {Q^c + Q^) cos9/N^cN^, J 



(11) 



where 



Q 



u 



Qd — 

Qu" — ~1 ~ 5^ 

Qd^ = l + fe 

= -|e 

= l + |e 

Qi/c = — 1 



(12) 



The factors Nf = ^1 + T^\Qf\'^ come from doing the algebra exactly rather 
than evaluating the lowest order graph in Fig. 1. [The linear combination of 
16's which has a Dirac mass with 16 and is thus superheavy is clearly seen 



So the superheavy u, for example, is {uiq + TQ^Uiq^)/ ^1 + T'^\Qu\'^, etc. 
Finding the orthogonal (light) linear combinations and writing the term 
QilSilSlOn which appears in Eq. (1) in terms of them, one obtains the 
exact expressions in Eqs. (8) - (11).] 

The striking feature of Eqs. (8) - (11) is that the mass matrices are rank 2. 
This is a consequence of 'factorization';^'"^ that is, that the mass matrices do 
not come from a Yukawa coupling coefficient that is a matrix in family space 



from Eq. (1) and the definition of T to be (16 + TQi63l63)/^l + T^jQ 
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but from a product of Yukawa coefficients that arc vectors in family space. 
The rank 2 comes directly from the fact that X]j(ai 16^)16 10 h involves just 
two distinct linear combinations of the light generations, namely ciilSi and 
16. In other words, the fact that there are two heavy generations is a direct 
consequence of the fact that Yukawa terms are bilinear in matter fields. (It is 
amusing to note that in this context if there were less than three generations 
then there would be no light family such as makes up the ordinary matter of 
our world. Perhaps this is a partial answer to the famous question of Rabi. 
Without the // and the r one would not have the electron!) 

A second striking feature is that these equations provide an explanation 
of the puzzling fact that the minimal SU{5) prediction^ = m° works so 
well while the corresponding predictions for the lighter generations fail badly. 
(Here and throughout, the superscript ° refers to a quantity evaluated at the 
GUT scale.) m° = m° is a consequence of the fact that D^s = L33, which 
follows from the relation Qd^ + Qd = Qi+ + Qi-- This relation, in turn, is 
implied by the fact that both d and / get mass from the same Higgs field, 
H' , so that Qd'^ -\- Qd + Qh' — Qi+ + Qi- + Qw — 0. On the other hand, 

- ~Li^ - -Di^' SO tnat m^/m, - -^^^ - -q^, wtiicn 

is a group-theoretical factor of order unity but not equal to unity in general. 

Another way to understand the fact that = while 7^ m° is 
from the group-theoretical structure of the term 16^16^10/745^:/. {Cf. Eq. 
(4).) The IOh x 45// can be contracted into a 10, 120, or 320 of 50(10). 
320 is not contained in 16 x 16 and so does not couple, li i — j — 3 then 
the contraction 120, which must couple antisymmetrically in flavor, is also 
not allowed. Hence only the contraction into a 10 couples to I63I63, and 
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therefore the minimal SU{5) relation D33 = L33 holds. But for = (2, 3) 
or (3, 2) the antisymmetric contraction into a 120 also couples, and thus the 
minimal SU (5) relation does not hold for the second generation. 

Three more facts are explained if the assumption is made that |e| <S 1, 
ie. that (45^) points approximately in the I3R direction. First, since 



^Q^^Q~ — 3 Y^j^ , the Georgi-Jarlskog relation^ m° = 3m° is explained. 
Second, since Qu, Qd, and Qi- are all proportional to e {ul, dL, and are 
all singlets under SU{2)ji), the matrices U, D, and L all become rank 1 in 
the limit e — >^ 0. Thus the hierarchy between the second and third generation 
masses is explained. It is remarkable that a relation among generations is 
related to a direction in the gauge group space. And, third, the angle Vcb is of 
order e. {Vg = (it) - = tan^ - ^] - tan^(g, + Q„) = 

|etan^.) This explains the smallness of Vch compared to ^Jm^/m^, which is 
a puzzle in Fritzschian models.^ In fact, since wPg/m^ = |esin^ 9, one has for 
e ~ 1 that - 2 (^^^ , which is a true relation. 

A further consequence of the structure of Eq. (4) is that the Higgsino- 
mediated proton-decay amplitude is proportional to^ e. Thus the smallness 
of e helps^ solve the problem of excessively rapid proton decay which tends 
to afflict supersymmetric grand unified models,^ especially those with large 
tan/9. (It should be noted that many SUSY GUT models based on 5*0(10) 
give tan/? = m°/m°.) The branching ratio Br[{p K^V^)/{p 7r"'"I7^)] is 
a calculable quantity in the modeP, but it differs from that of minimal SUSY 
SU{5). Similalry, the charged lepton mode p — > K^/j,'^ which could become 
significant in the large tan/? scenario has a different branching ratio than in 
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SU{5). These distinctions could serve as testing grounds for the high scale 
flavor structure. 

In summary up to this point, the simple structure of VFspinor together with 
the assumption that Q ~ I^r has explained or contributed to explaining six 
facts. And the single group-theoretical assumption about the direction of Q 
has played a role in four of these explanations, a striking economy. It should 
be noted that the direction I^r is not an arbitrary one but a point of higher 
symmetry, so that certain simple superpotentials can give an adjoint Higgs 
a VEV in that direction. For example, a superpotential with the terms 
W = + + SA^ + A^, where A is a 45 and S is a 54 has solutions 
for A in the I-i,r and B — L directions. And a superpotential with terms 
W ^ A^ + {A^y + A'^ has solutions in the I3R, B - L and X directions'^ 
{X is the S'C/(5)-singlet generator of SO{10).) It is also interesting that the 
Dimopoulos-Wilczek mechanism,'^ which appears to be the only viable way 
to make the gauge hierarchy natural in 5*0(10), requires there to exist an 
adjoint VEV in the B — L direction. 

At this point, more needs to be said about the factors Nf = ^Jl + T^IQ/p- 
Note that for T small, which corresponds to small mixing between the 16j 
and the 16, these factors become very close to unity and can be ignored. [T 
cannot be too small, however, since m° = aTvcosO, and therefore a > 1/T. 
For a ^ 1 perturbation theory would break down above the GUT scale.] 
Even with T of order unity the Nf become simple in the limit | e | <^ 1 , where 
N^^^^i- = 1 + 0(e2), and N^c^a-,i+ = v^TTT^ + 0(e). Effectively, then, there 
is for small e only a single parameter, N = Vl + T'^ ~ 1, introduced by these 
mixing factors. This parameter, iV, plays no significant role. For without it 
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there would be three parameters (e, 6*, and tan/5) to fit the five mass ratios 
and one mixing angle of the second and third generations. That means there 
would be three quantitative "predictions", which are (for the small e case) 
m° = m°, m° ~ 3m°, and m°/m° = m°/m° (about which more later). But 
one easily sees that the factors of appearing in the (32) and (33) elements 
of all the mass matrices approximately cancel out in all these relations and 
thus leave them unaffected. As for the two qualitative "predictions" , namely 
that Vcb and the ratios of second to third generation masses are of order e 
and hence small, these are clearly unaffected by the presence of N, which is 
of order unity. 

The simple structure of W^spinor (Eq. (1)), as we have just seen, explains in 
a highly economical fashion most of the features of the masses and mixings of 
the two heavy generations, and indeed the very fact that there are two heavy 
generations. Two facts about the heavy generations are not explained. The 
first is that ^ m°, which is just equivalent in this model to the statement 
that tan/3 3> 1. An explanation of this must rely on an understanding of 
the Higgs sector. If we regard this model as a model of the Yukawa sector, 
this question is beyond its scope. The second unexplained fact is the failure 
of the "proportionahty" relation of SO{10), wP^/m^ — mPg/m'i,, which is not 
significantly broken by any effects discussed so far. Different mechanisms 
for breaking this bad relation lead to different versions of the model, see 
subsections (c) and (d) below. 

(b) The first generation: VFyector 

The masses and mixings of the first generation arise from the next layer 
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of the model, given by W^ector- As before, in the case where the mixing of 
the 16i with the 10 + 10' is small the contribution of H^vector to the light 
quark and lepton masses can be read off from a simple graph, shown in Fig. 
2. At first glance it might seem that many new parameters are introduced 
by adding VFvcctor- (In fact, nine: the six Yukawa couplings, q and c^, and 
three VEVs.) But the important point is that in the limit of small mixing 
Fig. 2 gives a flavor- antisymmetric contribution to Dij and Lij. (There 
is no contribution to Uij as 10 and 10' do not contain charge-(|) quarks.) 
The flavor- antisymmetry is obvious from Fig. 2, for under the interchange 
10 10' the Yukawa coupling 10 45^10' changes sign because the adjoint 
of 5*0(10) is an antisymmetric tensor, while the flavor indices i and j are 
interchanged. Moreover, if it is assumed that the 45h acquires an SU{5)- 
invariant VEV, the contributions to Dij and Lij will be equal (up to a sign 
since 5D — 6LF). Thus in the small mixing limit (we will return to the more 
general case later) there are really only three new parameters introduced 
which we will call Qj, {i ^ j)i where 

^^^M(45^)(45^)(5(10^))- ^''^ 
It should be noted that if the Cij are all comparable then C23 is small compared 
to the other contributions to the terms in which it appears {Diz-, -D32, L23, 
and L32) and can therefore be neglected. 

The mass matrices then take the form (putting in also the values of the 

from Eq. (12)) 
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U^aTv 



( 









(14) 



V -(1 + |e) sinO/N^c -(1 + |e) cose/N^cN^ 



D = aTv' 



-ci2 -cn/Nd \ 

ci2 (lesmO -C23)/Nd 

V ci3/iVd^ ((1 + |e) sin ^ + C23)/A^dc (1 + |e) cos e/NacNa J 



(15) 



aTv' 



( ci2 CiaM- \ 

-ci2 (-fesm^ + C23)M- 

V -C13M+ ((1 + |e) sine - C23)M+ (1 + |e) cose/Ni+Ni- ) 



(16) 




N = aTv\ Q 





-|e sin 6'/iV^ 



(17) 



-sinO/N^c -(l + ie)cose/iV^ciV^, 



Several more features of the quark and lepton mass spectrum are ex- 
plained by this form. First, that Uij is still rank 2 corresponds to the fact 
that m^/m^ (~ 10~^) is much smaller than wPjm^ (~ 10~^) and w^l/m^ 
(= 0.3 X 10~^). Second, the antisymmetry of the contributions to the first 
row and column of implies that one has effectively the form ^ ^ ^ 
for the (12) block of that matrix (after diagonalizing the (23) block). This 
is precisely the well-known form^^ that leads to the famous relation^^'-^^ 



tan 6'c = y'm^/m^. (There is no y'm^/m^ contribution. Some other ef- 
fect will generate a non-zero and very small m„, but there is no reason for 
this other effect to have a form that makes C/n <^ C/12, C/2i- Rather, if one 
supposes that C/n, and U21 are all of the same order then the contribu- 
tion from diagonalizing U to tan 6'c will be 0(m°/m°) ~ 0.005 instead of the 
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usual ym^/mj? = 0.07. Thus this "prediction" works better here than in the 
usual scenarios.) 

A third consequence of the forms given in Eqs. (15) and (16) is that 

detD = detL. (18) 

Note that, remarkably, this is true for any values of ^, e, and Cy , if T is small 
enough that the factors Nf can be taken to be one. And for any T, this 
equality holds for small e. Prom this follows the well-known Georgi-Jarlskog 
relation^ ml/wP^ = (m°/m°)-^ = i. 

Finally, if all the c.^ arc comparable then Vub ~ VusVcb, which is the true 
order of magnitude statement which underlies the Wolfenstein parameterization. 
(For Kb ~ Dsi/D^s ~ aTv'{ci3/m^^), ~ aTv'{ci2/m'',), and Va ~ mO/mO, 
the last being both empirically true and a consequence in the model of ^ ~ 1 
as noted earlier.) 

(c) The proportionality relation: The short model 

The structure described so far, given in Eqs. (1) and (2), gives a sat- 
isfactory account of all the features of the quark and lepton spectrum (the 
threefold hierarchy, the mixing angles, and the mass ratios) with one glaring 
exception. The ratio m°/m^ is seemingly predicted to be equal to m°/m°, 
whereas empirically it is about one-fifth of that. The most obvious expla- 
nation of this anomaly would be that another additive contribution to Uij 
exists which happens to approximately cancel C/23- It is plausible that this 
might happen without notably changing any of the other successful features 
of the model since U23 is a small element and also does not affect the mixing 
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angles. This is an attrctive possibility. But it suffers from the apparent dif- 
ficulty that the new contribution to Uij must leave it very nearly rank 2 in 
order not to produce a value of m„ that is too large. We have found a way 
to do this, which is discussed in Section 6(b). This idea, however has certain 
drawbacks, discussed in Section 6(b), that make it less attractive than the 
ideas we shall now present. However, we cannot discount the possibility that 
a more elegant solution along these lines may exist. 

A remarkable fact about the model described so far, containing only the 
Yukawa couplings in Eqs. (1) and (2), is that it already contains, without any 
further additions, a multiplicative correction to the mass matrices that can 
break the bad proportionality relation while leaving the other good relations 
largely intact. 

We have emphasized that the forms of the mass matrices given in Eqs. 
(14) -(17) are valid for small mixing. The exact expressions for the D, L, 
and N matrices including all mixing effects are 



D = aTv'{I+Adc)-^ 



( -ci2 -cu/Nd 

Ci2 (|esin^-C23)/iVd 

V cis/iVdc ((1 + |e) sin 9 + C2i)/Nd^ (1 + |e) cos O/Nd^Nd 

(19) 



) 




(20) 
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/ 

N ^aTv 






sin 





-|esin^/iV^ 
(l + |e) cos 9 /N^cN, 



) 



(21) 



where 



{Adc)ij — C*Cj + C*'Cj, 



(22) 



and where 



(Ci, C2, C3) 



(1(16h)) 
d{45H) 



{ci,C2,C3/Ndc), 



(23) 



with a similar expression for C- in terms of the c-. The expression for A^- 
is of the same form with A*";- instead of Ndc, and Ajy = Ai- by SU{2)l. See 
Appendix A for the derivation. The expression for Uij given in Eq. (14) is 
exact as it stands. 

The matrices Af characterize the mixing between the IGj and the 10 + 
10', just as TQie^ characterizes the amount of mixing between the 16j and 
the 16. The factors (/ + Af )~^ are completely analogous to the factors of 
(1 + r2|Q/|2)-2 = Nj\ For small mixing these factors are close to the 
identity and can be neglected. Even if the Af are not small, most of the 
quantitative and qualitative successes of the model are only slightly affected. 
Obviously the fact that m„ ^ and the threefold hierarchy are among these. 
Moreover, if Nj = 1 then for arbitrarily large Af the relation detD = det L 
continues to hold, since then A^c = Af- . 

Consider a case in which 




(24) 
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where 6 < 1. From Eqs. (19) and (20) it can be seen that the main effect of 
this factor is to suppress the r and b masses by a factor of approximately S. 
This suppression can be understood simply and intuitively in the following 
way. 6 arises, as can be seen from Eqs. (22) and (23), from a large value 
of C3 or C3 compared to (i(45i/)/ (1(16^/)). From Eq. (2) one sees that there 
is a superheavy mass term 5(10) [(^(45^)) -5(10') + (03(1(16^))) -5(163)]. 
Large C3, therefore, corresponds to the superheavy linear combination being 
approximately 5(163) and the orthogonal, light linear combination that is 
the third generation being approximately 5(10'). That means that the term 
acos^^ 5(163)10(16) (5(16//)) in Eq. (1) gives a contribution mostly to the 
superheavy 5 rather than to the third generation 5 which contains b'^ and r~. 
Thus the masses of b and r are suppressed by a factor of approximately S. 

It is easily seen by substituting Eq. (24) into Eqs. (19) and (20) that the 
masses of b and r are multiplied approximately by 5 while the masses of s 
and fi are not much affected. Thus the ratio m^/m^ is enhanced relative to 
m^/m^ by a factor of 6~^. Hence 6 should be about |. At the same time, the 
near cancellation between the (23) mixing angles in the up and down sectors 
that is responsible for the smallness of Vcb is not disturbed, since the ratio 
-D32/-D33 is left unchanged. 

What this shows is that the factors (J + A)^^ in Eqs. (19) and (20) 
allow a fit to be made to {m!^/m!^)/{m^/ml) without greatly disturbing 
the qualitative and quantitative successes of the model. However, there 
are certain prices to be paid for this method of breaking the proportion- 
ality relation. First, 9 must be somewhat small. Since Vcb — fetan^ and 
m^Jml ^ ^esiir^ 9/S, it follows that 9 ^ 25^^^ - 5 ~ i. There is noth- 
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ing in principle wrong with this, except that as 9 is the angle between two 
supposedly unrelated vectors, a and 6, it might have been expected to be 
closer to unity. Second, the smallness of {m^^/mf) /{wPg/m^) is not so much 
explained as fit, and, in particular, a special choice of the form of (/ + A(jc)~5 
has to be made. {Cf. Eq. (24).) This choice amounts to the statement that 
the vector q (or c'j) points nearly in the '3' direction. There is then a kind 
of preferred direction in family space, in spite of our having eschewed fam- 
ily symmetry. (Though, on the other hand, the 'small numbers' involved in 
this fortuitous alignment are only of order |, whereas the intergenerational 
hierarchies being explained involve ratios of 10^^ to 10^^.) A third cost is 
that the exactness of the relation m° = m° is lost. For example, the simple 
form given in Eq. (24) would lead (for small e and C23) to = S{aTv') and 
m° = V 5'^ cos^ 6 + sin^ 9{aTv'). Since both 5 and 9 are of order |, m!^/m'^ 
can deviate from unity by as much as 40%. As we shall see in Section 3, the 
form of Adc can be chosen to fit m^/m^. But now it can only be claimed as 
a "prediction" that the masses of b and r are equal in a rough sense. Finally, 
because the directions of Cj and are not arbitrary but must be chosen to 
give Arfc the required form, it turns out that there is a hierarchy among the 
Cij. In particular, C23 is somewhat larger than C12 and C13, and therefore, 
while still small compared to the terms in which it enters, it is not so small 
as to be negligible. In fact, the Georgi-Jarlskog factors of 3 are significantly 
affected. (See Section 3.) 

(d) The long version of the model 

While all of the costs of fitting {m^^/m^)/{mym^) that had to be paid in 
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the short version of the model are relatively minor, they do somewhat take 
away from the cleanness of the model and its explanatory power. There is, 
however, another very beautiful way to break the proportionality relation. 
This involves adding the terms Wadoint given in Eq. (3) to the superpotential. 
To some extent this reduces the economy of the model, though as emphasized 
earlier the overall structure of VFadjoint is parallel to that of PVspinor and W^vector- 
Moreover, by adding these terms one is enabled not merely to fit but to 
explain the smallness of {my m^) / {mf^ / m'^) in a rather elegant way, as will 
now be explained. Finally, this explanation is achieved in such a way that 
the other positive features of the model are virtually unaffected. None of the 
costs that have to be paid in the short version have to be paid here. This 
then is a much cleaner version of the model. 

If one assumes that (5(16ij)) = (this will be seen to be desirable on 
other grounds in Section 5) there is no additive contribution to the mass 
matrices (analogous to the q^) coming from VTadjoint- W^adjoint does introduce, 
however, new multiplicative corrections of the form (I+Af)~^, where / = u'^, 
u, d, and Z"*", that reflect the mixing of the 10(16j) with the 10(45) and 
10(45'). This mixing comes from the following terms contained in Wadjoint^ 

T0(45) [(/(45^^)) • 10(45') + E^(e^(l(T6^))) • 10(16,)" 
_ _ (25) 

+ 10(45') [(/(45^)) ■ 10(45) + E.(e^(l(16//))) ■ 10(16,) . 

The effect of this is to introduce factors of the form (7 + A/)~^ into the mass 
matrices U, D, and L wherever they multiply an 5'C/(5) 10 of fermions. For 
example, the matrix Uij now takes the form 
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/ 



U = aTv{I+Auc)-^ 








V -(1 + |e) sin^/iV„c -(1 + fe) cos 6 /Nu^Nu 



The matrices A„c, A„, A^, and A;+ all have the form 



(/+A, 
(26) 



(A;),, = e:e,+ Er^, 



(27) 



where 



{El, E2, E3) 



(1(16^.)) 



(ei,e2,e3/Arj), 



(28) 



/f2Q/(45) 

and similarly for the definition of E'-. (The expressions for A^c, A;-, and A^, 
have already been given in Eqs. (22) and (23).) 

A point of crucial importance is that the expression for Aj, / = u'^, u, d, 
and given in Eqs. (27) and (28) has a factor of Q/(45) in the denominator. 
This simply comes from the VEV of the 45 h that appears in Eq. (25). Note 
that this adjoint is the same one that appears in H^spinor (rather than the 45h 
that appears in Wyector) and has its VEV in the direction Q. But it can be 
seen from Eq. (3) that the generator Q in VFadjoint is not acting on fermions 
that are in a 16, as in VTspinoD but on fermions that are in a 45. What is the 
difference? Using the relation /3R = —j^X + | (^)' rewrite Eq. (7) 

as 

Q=(-i)x+(|(l + 6))f (29) 

The SU (5)-singlet generator X (as we have normalized it) takes the value 1 
on the representation 10(16), but -4 on the representation 10(45). Thus 



<3/(45) = <3/(16) + 1- 



(30) 



18 



Thus the charges that enter into Eq. (28) are 





(31) 



Note the important fact that for |e| <^ 1 the only field in the 45 that has 
a small charge is the w"^. The remarkable consequence of this is that the 
elements of A„c are of order ^ (see Eqs. (25), (26), and (29)), while all 
the other Af have elements that are not enhanced in this way. In this long 
version of the model, then, the assumption can be made that all the A/ are 
less than unity so that the factors (/ + A/)~^ are unimportant, except for 
the matrix A„c, which has large elements of order ^. 

Because of the foregoing only the matrix U gets substantially affected by 
the presence of H^adjoint- Most of the successful relations of the model are 
then essentially unaffected, since they follow from the forms of D and L. 
The quantities that depend strongly on the form of U are m^, mc/mt, and 
the mixing angles. But, as the factor (/ + A^c)"^ does not alter the fact that 
U is rank 2, the smallness of m„ is still explained. And because the factor 
(/ + Auc)~5 involves a transformation of the right-handed quarks, the mixing 
angles are not strongly affected. One might expect that the cancellation 
between the contributions from the up and down sectors that makes Vcb 
small would be seriously disrupted if (/ + A„c)^5 were very different from 
the identity matrix. It turns out, however, that Vcb remains of order e. In fact 



Prom these considerations it is seen that the only prediction of the model 
that is substantially altered by the addition of the term VFadjoint is the pro- 



(see Appendix B for the derivation) 1/° = 



lesin^cos^ l + ftan^j|^ 

5 2 {ExE')3 
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portionality relation for m^/m^. At first glance one might expect that the 



m'^/mf, for arbitrary values of the parameters in A^c But, remarkably, this 
proves not to be the case. In Appendix B we give an explicit proof which 
shows that if the elements of A^c are large then m^/m^ will be suppressed 
except for special directions in parameter space. (Interestingly, this is a con- 
sequence of the fact that A„c is a rank 2 matrix, as can be seen from Eq. 
(27). The ratio {m^/m^)/{rnPg/m1) in this case is given by 



from which is is clear that {m^^/m^) is suppressed relative to {m^g/m^) for 
\Ei\ > 1 (i.e., e < 1, Cf. Eq. (28)). See Appendix B for a derivation. If it 
were an arbitrary matrix of rank one or rank three, then m^/m^ would with 
equal likelihood be suppressed or enhanced. This also is shown in Appendix 
B. These statements have also been checked by numerical tests.) Of course, 
if the elements of A^c are small, m°/m° will be only slightly affected one way 
or the other. In particular, since the elements of A^c are of order ^ the ratio 
m°/m° is suppressed by a factor of order e as shown above (see Appendix 
B), which is what is needed to agree with experiment. 

The surprising conclusion is that the addition of the terms VFadjoint to 
the superpotential for generic values of the parameters suppresses the ratio 
m°/m° by order e, while having only a minor effect otherwise. The violation 
of the proportionality relation is therefore not merely fit in some arbitrary 
way, but explained in a group-theoretical way. In this version of the model, 
the relation |e| <^ 1 plays a role in explaining no less than five facts! Such 



factor (/ -I- A„c) 2 would be just as likely to increase as to decrease the ratio 
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economy of explanation is not something that can be contrived. 



3 Numerical Fits 

(a) Long version 

In this Section we shall present the details of our numerical fits. Let 
us first focus on the long version of the model. In this version, we can 
consistently set the multiphcative matrices A„, A^, A^c, A/- , A,^ and A/+ (Cf. 
Eqs. (19)-(21)) all to zero. The elements of A„c are enhanced by a factor 
l/|e|^ and therefore A„c cannot be ignored (Cf. Eq. (26)). In this limit we 
have the following approximate analytic expressions for the mass ratios and 
mixing angles: 



m. 







aTv \{E X E')3 
N \Ex E'\ 




fesm^9\E x E'l^Ef + E'-^ 



{E X E%\^ 



~ 



(32) 




Ci2-^^(ci2COS^ — Ciasin^) 



(33) 







~ 3m°, 




(34) 



V 
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V^^ ~ -ecos^sin^ 



us — 



iV X #)2 

2"(Ex#)3 



cos6 — — sin^ 

Ci2 



-1/2 



(35) 



Note that in this version, the angle 9 is of order one, as a result, there 
is a significant correction to the expression tan^c — However, the 

ratio C13/C12 is nearly equal to 0.5 in order to fit V^b, and with this value the 
correction factor in Eq. (25) is small for a wide range of the angle 9 as long 
as Ci3 and C12 have a relative negative sign. (For example, if C13/C12 = —1/2, 
the correction factor in Eq. (25) is 1.04 corresponding to ^ = 60°.) This fact 
is borne out in our numerical fits. 

A good fit follows by choosing 



T 

C12 



e = 0.15i, e= (1,1,0), = (1,-1,0) 

4fn 



1, sin^ = cos^. 



= 1.2 



5 (45h) 

0.0035, ci3 = -0.002, C23 = 0.002 



(36) 



The resulting mass eigenvalues are 

(m°, m°, m°) = [0, 3.64 x 10'^, Q.7Q8]aTv 

m°, ml) = [8.54 x 10'^ 1.59 x 10-^ 0.709]aTv' 
{ml, mj, m°) = [2.98 x 10-^ 4.53 x 10"^ 0.707]aTv 

The quark mixing angles have the values 

K° = 0.213, K° = 0.002, V;° - 0.030 



(37) 



(38) 



Since the input parameters are all complex in general, there is room for 
sufficient KM type CP violation. In the specific example given above, this 
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CP violation is somewhat small, but for other choices of input phases large 
enough CP violation can be obtained, 
(b) Short version 

Here the fit is somewhat nontrivial since the same set of vectors q and c[ 
of Eq. (2) should generate the antisymmetric contribution of Eq. (13) as 
well as the multiplicative factor A^c = A^- given in Eq. (19)-(20) needed to 
correct the proportionality relation. As already noted, this can be achieved 
by choosing an approximate form of A^c as given in Eq. (24). However, 
and m° are now only approximately equal for generic values of the model 
parameters, to within 40% or so. We shall choose the parameters so that 
m° ~ mP^ remains good to within about 10%. The angle 9 in this case is of 
order 1/5, so that the correction to the relation tan^c — ^/^^^^/^^^^ is small. 
We will also need T < 1 so that ~ 1. 

The following choice of parameters gives a good fit: 

c = (1,0.7,7), (? = (0.75,1.5,4.5) X 10"^ 

= {cic'j-c'^Cj), sine = 0.22, e = 0A5i, T = 0.25 (39) 

The masses and mixing angles with this choice are 

(m°,m°,m°) = [0, 4.5 x 10-^ 1.0]aT^; 

{m% m°, ml) = [1.7 x 10'^, 4.2 x 10'^, 0.22]aTv' 

(m°,mj,m°) = [5.2 x 10-^ 1.2 x 10"^ 0.25]arT;' (40) 

C = 0.213, V^, = 0.0025, Kb = 0.0295 (41) 
Note that these numbers correspond to the ratios m^^/mP^ = 2.82, my m'^ = 
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3.26, m°/m° = 1/225, m°/mg = 1/52, (m°/m°)/(m°/mg) = 1/4.3, all of 
which are in good agreement with data. 

4 Predictions 

(a) tan/3 

If one neglects all the factors of {I+Af)~^, then v/v' = m^/m1 = m^^/m^g. 
Of course, it is the failure of the latter equality that requires one to assume 
that not all of the mixing matrices are trivial. In the short version of the 
model, the mass of the h quark is suppressed by such effects while the masses 
of s, c, and t are little affected. Thus v/v' = m^/m^ still holds but m^/m1 
is larger. If one assumes that (5(16if)) <^ {5{10h)) (see the discussion in 
Section 5(a) for why this is sensible) then tan/3 = v/v' and there follows 
the unusual prediction that tan/5° = m^/m^^. After renormalization group 
corrections, tan/3 at weak scale is w 18 (corresponding to rrit — 175 GeV .) 

In the long version of the model, the masses of h and s are little affected 
by the mixing matrices, but both and are suppressed. (See Appendix 
B.) However, m° is only suppressed by a factor of order unity, while is 
suppressed by a factor of order e. Thus (again assuming tan/3 = v/v') one 
has the prediction tan/3 ~ m^/ml. (This relation is unaffected by renormal- 
ization group running.) 

(b) Neutrino mixing angles 

The predictions for neutrinos will be discussed first in the long version 
of the model where they arc simpler. In the long version it is assumed that 
all the A/ are somewhat less than unity so that the factors (/ + A/)~^ are 
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close to unity and negligible, except for A^jc, which is enhanced by the group- 
theoretical effect discussed above. But obviously the factor (/ + A„c)~5 does 
not affect the neutrino mixing angles at all. 

At first glance it might be thought that there could be no predictions 
for neutrinos in cither version of this model, since the neutrino mass matrix, 
My, depends not only on the Dirac neutrino mass matrix, A^, which is known 
(after fitting the charged fermion masses and mixing angles), but also on the 
superheavy, Majorana mass matrix, Mr, which is unrelated directly to the 
others and not known or predicted in this model. 



However, ignoring those effects which produce m„ 7^ 0, the matrix N , like 
[/, has rank 2, with vanishing first row and column in the basis in which we 
have worked. (See Eq. (11).) It is then clear that (in this hmit) for any 
form of M^^, the matrix M^ also has vanishing first row and column. Thus, 
in diagonalizing M^ no rotation in the 1-2 or 1-3 planes is necessary, and 
the mixing angles Ve^ and Ver come entirely from the diagonalization of L, 
which is a known matrix, just as \4s and Vuh come from the diagonalization 
of D. Since the form of L is similar to that of D, there are relations between 
the lepton and quark mixing angles. In particular, since L33 = D33 and 



My = N^M^^N. 



(42) 



= — £>3i, one has ^^3'' ^ _ 



(Vkm)i3, or 



— —Kb- 



(43) 



And for the same reason that \4s — 





(44) 



25 



It is likely that whatever effect makes U a rank 3 matrix also makes N 
rank 3. Just as these effects are likely to correct tan^c by 0(m°/m°), one 
would expect that the corrections to and Ver would be roughly of order 
m1/m^ ~ 0.005 and m^/m^ ^ 10~^ respectively. But it is possible also that 
the effects that produce m„ act differently on the neutrinos. 

The mixing angle Vfj_r cannot be predicted in the absence of any infor- 
mation about Mr. However, in certain limits interesting predictions arise. 
Simply multiplying out Eq. (42) one finds that 

(M,)33 = {M],'Ml + lefcos'e/N' 

+ 2(M^i)23fe(l + |e)sin^cos^/7V 



(M,)23 = {M,h2 = iMj,')s3{l + l€)sm9cose/N' 



(45) 



+ {M^')23hsm^e/N, 



(M,)22 = iM^')3ssm^9/N\ 
And, of course, (Mj,)ij = (Mj,)ji = 0. From these equations it follows that 



(M^i)33 > (M^i)23, e(M^i) 
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l + 0(g|;gtan^)+0(6: 



(46) 



Another interesting case is 

(M^^)23 > (etan^)-i(M^^)33, etan^(M^^)22 



(47) 



m{v^)/m{vr) = i- 

In the short version of the model the neutrino predictions are complicated 
by the presence of the factors (/+ A;-)~^ — (7-|- A,^)""^ in Eqs. (20) and (21). 
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These factors affect the mixing of the left-handed charged leptons. In the 
absence of these factors one would have the same predictions for neutrinos 
as in the long version. 

The matrix A;- is given by the expressions in Eqs. (22) and (23) with Nd,c 
replaced by A*";-. A satisfactory suppression of (and thus of (m^/m^) / (m^/m^)) 
would be achieved if ci, C2, c'j^, and were very small compared to one, 
while C3 and/or C3 were larger than one. For simplicity of discussion let us 
henceforth ignore c\ and just assume that C3 > 1. Then (Ai-)33 = \c2,f > 
1 and all other elements of (A^-) are small, which gives (/ + A;-)~5 ^ 



diag(l, 1,1/ Jl + This has the effect of multiplying L33 by 



and thus the rotation in the (yU^, r^)-plane required to diagonalize L is not 
approximately tan^ but ^^1 + jcal^tan^ = ^0/^0 tan 6, which numerically 
(see Section 3) is of order unity. The short version, then, typically predicts 
that V^T is large. 

A second effect of the matrix (/ + A;-)^5 arises from the fact that Ci 
and C2 are not in general exactly zero. In fact, as was seen in Section 3, 
one can fit the relation = m° by choosing Re(c2/m3) = — |tan^. Now, 
the presence of a non- vanishing ci, which would be natural to expect, leads 
to non-vanishing (12) and (13) elements of (/ -|- A;-)~^. These, in turn, 
contribute to and Ver- In fact, Ver can be quite large. Because of the 
presence of the unknown parameter Ci, there are not independent predictions 
for Ve/j, and Ver, but one prediction for Ve^ in terms of Ver- For Ver small, the 
prediction for Ve^ just goes over to that for the long version of the model. 
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More generally 

Km = V^°/mO + 0(e). (48) 

5 Technical Issues 

(a) The Higgs sector: doublet-triplet splitting 

In 5'O(10) SUSY GUTS the simplest way to naturally achieve the doublet- 
triplet splitting of the Higgs fields that is required for the gauge hierarchy 
is by the Dimopoulos-Wilczek mechanism. The essential idea involves an 
adjoint of Higgs (^41) whose VEV is in the B — L direction and a pair of 
fundamental Higgs multiplets (Ti and T2). Consider the following Higgs 
superpotential. 

1^2/3 = ATiAiTs + M^iT^f. (49) 

Since the Tj are in 10 representations, each contains a 5 + 5 of SU (5) that 
we will denote 5(Ti) and 5(7]). Then one has the following mass matrix for 
those fields 

If {Aij = a - {B — L) + b- (/sr), then, in an obvious notation, one has for the 
masses of the 5't/(2)L-doublets and color-triplets contained in Tj 

■ Xb \ ( 2(Ti) \ 



W^mass = (2(Ti), 2(7^2)) 



+ (3(ro,3(r2)) 



-Xb M2 )\ 2(T2) 
Aa A / 3(ri) 



(51) 



-Xa M2 J \ 8(^2) J ■ 
If a ~ Mqut ~ M2, while 6 = 0, there is a massless pair of doublets, 
2(Ti) + 2(Ti), that play the role oi H + H' in the super symmetric standard 
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model (SSM), while the other pair of doublets and all of the triplets be- 
come superheavy. Moreover, if one assumes that T2 does not couple to light 
quarks and leptons the proton-decay amplitude coming from the exchange of 
the color-triplet Higgsinos, 3(Ti) -|- 3(Ti), is proportional to M2/{Xa)^. For 
the theory to be perturbative A cannot be large compared to unity, and so 
Xa ^ Mgut- However, M2 can be somewhat smaller than Mgut, and this 
provides a means by which the proton lifetime can be made consistent with 
experiment.^ (In the minimal SU{5) SUSY GUT it is well-known that the 
proton lifetime is only marginally consistent, if at all, with experiment.^) 

In general, b will not stay exactly zero if non-renormalizable terms (in- 
duced by gravity, for example) are taken into account, and, in fact, it is 
not a completely trivial matter to ensure that nonrenormalizable terms do 
not destabilize it. How small must b remain to preserve the gauge hierar- 
chy? As noted, the proton-decay amplitude is proportional to M2/(Aa)^, 
which must therefore be < (IOMqut) The mass of the light Higgs dou- 
blets, on the other hand, is given by (A6)^/M2, and therefore can be written 
niH ^ (10 Mgut). thus, it must be that b/a < 10"^ < (Mgut/Mpi)^. 
That this can naturally be achieved is shown in Ref. 16. 

This simple picture of doublet-triplet splitting must be modified in the 
context of the present model, since in addition to the IOh appearing in 
W^spinor, which is to be identified with Ti, there is the 16h appearing in 
M^vector, which also must break SU{2)l x U{1)y and must as a consequence 
partially contain a hght 2. To be more exact, there is just a single hght 2 that 
gets a Weak-scale VEV and is a hnear combination of doublets in 10// and 
16h- This raises two issues: (1) Can the Dimopoulos-Wilczek mechanism 
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for doublet-triplet splitting still work?^^, and (2) does the mechanism for 
suppressing Higgsino-mediated proton-decay still operate? 

In the long version of the model there is also a, 16h that is required 
not to get an SU{2)l x C/(l)y-breaking VEV. So issue (3) is whether this 
is natural. Finally, it is desirable that (5(16H))/(5(10i/)) (which we shall 
define to be tan 7) be small for two reasons. First, it would explain (see 
Eq. (13)) why the are small, and, second, it would mean that the light 
Higgs, H', of the SSM would be almost purely in the 10//, so that the ratio 
(5(10/f))/(5(10//)), which is predicted in the model, is just the empirically 
measurable parameter tan/?. Issue (4) is whether the smallness of tan 7 can 
be simply and naturally achieved. 

Let us denote IOh, 16//, and 16// by Ti, C, and C. A satisfactory 
generahzation^''' of Eq. (49) is 

1^2/3 = ATi^iT2 + M2 (T2)2 + pT^CC + McCC. (52) 

which gives 

/ X{A,) p(C) \ I 5(Ti) \ 
Vr^ass = (5(Ti),5(T2),5(C)) -A(Ai) M2 5(T2) . (53) 

V Mc I \ 5(C) / 

If 6 = (that is, if {Ai) = a{B — L)), then the matrix for the doublets has 

one massless eigenvalue. The massless 2 is purely in Ti, while the massless 

2 is a hnear combination of 2(Ti) and 2(C): 

2iight = cos72(Ti) + sin 72(C), (54) 

with 

tan7 = -p(C)/Mc. (55) 
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Since the orthogonal, superheavy doublet, 2heavy = — sin72(Ti) + cos72(C), 
must have vanishing VEV, it follows that (2(C)) /{2{Ti)) = (5(16i/))/(5(10j^)) 
= tan 7. Thus by making the ratio p{C)/Mc small one ensures that issue 
(4) raised above is satisfactorily resolved. One possibihty, discussed later, 
is that the term pTiCC arises from a higher- dimension operator, so that 
p~ 0(Mgut/Mpi). 

It is necessary, as before, that h/a < 10~^ to preserve the gauge hierarchy 
(issue (1)), but it is no longer sufficient. It is also necessary that the lower 
left entry in Eq. (53 ), call it x, which connects 5(C) to 5(Ti) be extremely 
small. To be precise, it must be that xp{C) / AIq ^ iny/i or, in other words, 
X ^ cot ^iiTLw- It should be noted that this would also automatically ensure 
that (5(T6i/)) ^ 0, which resolves issue (3). For (5(T6j^))/(5(10i/)) = 
(2(C))/(2(ri)) = x/Mc < mw/p{C) ~ mw/MovT. 

The remaining question (issue 2) is whether Higgsino-mediated proton 
decay can still be suppressed by making M2 somewhat smaller than Mqut- 
This is easily seen to be the case from the form of Eq. (53). If M2 is set 
to zero, then the colored Higgsino 3(Ti) (which is by assumption the only 3 
to couple to the light quarks and leptons) only has a mass connecting it to 
3(T2), which does not mix with the 3(Ti) and 3(C) that couple to the light 
quarks and leptons. Thus the Higgsino-mediated proton-decay amplitude is 
proportional to M2 as in the simpler case of Eq. (49). Note also that there 
is a group theoretical suppression of order e in the rate for proton decay. ^ 

In order for this simple scenario to work naturally, we have seen that 
several conditions must be satisfied, (a) As before, b/a must be < 10'"'''. (b) 
For a; to be sufficiently small the coefficient of any effective TiCC term must 
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be less than about mviz/AfcuT- And (c) the TiCC term should have a small 
(0(Mgut/-/W^pi)) coefficient. In the next subsection we will see that a realistic 
superpotential satisfying these criteria can be constructed. 

(b) The Higgs sector: the breaking of 50(10) 

The model of quarks and leptons requires the existence of adjoints with 
VEVs in the I^r and X directions. We will denote these by and ^3, 
respectively. The Dimopoulos-Wilczek mechanism requires the existence of 
an adjoint whose VEV is in the B — L direction, which we have been denoting 
Ai. (As above, we will denote the 16^^ and IQh by C and C, and the two 
lOif by Ti and T2.) A satisfactory form for the superpotential is 

W^Higgs = tr(Ai)7Afpi + Af^i(Ai)2 

+ tr(A2)7Mpi+_M^2(^2)' 

+ McCC + ^C{A^fC + MaM^Y (56) 

+ iY{AiA2A^)^' 

We have not written the dimensionlcss coefficients, which are assumed to be 
of order unity. This as a hybrid of the forms proposed in Ref. 16. 

The form W{A) = tr(A^)/Mpi + M(A)2 has as possible solutions A(xX, 
B-L, and I3R, with ~ ^/MpiM. Thus Mai and Ma2 in Eq. (56) must 
be of order M^urp/Mpi. 

The terms involving A3, C, and C can be shown to have a solution in 
which the VEVs of these fields are in SU (5)-singlet directions. The magni- 
tude of (As) is determined by the Fc and F-^ equations to be O(Mpi). This 
means that SO{10) is broken to SU{5) x U{l)x near the Planck scale. The 
fractional mass sphttings within SU (5) multiplets will then be 0(Mgut/A^pi) 
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and many of the threshold corrections to sin^ 9w will be suppressed by that 
small ratio. (See Ref. 16.) On the other hand, it is desirable that the VEV 
of C be of order Mqut because of the role it plays in Wadjoint- The scale of 
(C) and (C) is determined by the Fa3 equation to be Ma^/ McM^{). 

The term tr(AiA2A3), by linking the Ai, A2, and A->, sectors, ensures that 
there are no goldstone modes, while at the same time not destabilizing the 
VEVs of the Aj. (See the first paper of Ref. 16 for a discussion of this term.) 

The terms involving Tj have already been discussed in the last subsection. 
If {R) ~ Mgut, then tan 7 ~ Mgut/-/^pi as desired. 

(c) Discrete symmetries 

The essential core of the model of quark and lepton masses consists of 
the Yukawa terms given in Eqs. (1) - (3). In particular, the "root model" is 
defined by the terms in Eq. (1). The most important feature of that set of 
terms is that there is no direct gijlQilQjlQn coupling. This structure can be 
explained by a very simple Z2 parity under which the Higgs fields, lO^j = Ti 
and 45 H = ^2, and the families, 16j = Fj, are odd, while the extra real 
representations of matter fields, 16 = F and 16 = F, are even. 

The structure of the complete set of Yukawa terms given in Eqs. (1) 
- (3) can be ensured by a Z3 x Z2 symmetry, where Z^. is a matter parity 
under which "matter fields" are odd and Higgs fields are even, and Z3 is 
a symmetry which acts on the fields as shown in Table I. (The other extra 
real representations of matter are denoted 10 = T, 10' = T', 45 = A, and 
45' = A'.) 

This particular symmetry allows in addition to the terms in Eqs. (1) - 
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(3) a few extra harmless terms (namely FFTi, FTC, and FT'C). 

A realistic Higgs superpotential (that is, one which completely breaks 
SO{10), avoids goldstone modes, preserves the unification of couphngs, natu- 
rally achieves doublet-triplet splitting and sufficient suppression of Higgsino- 
mediated proton decay, and gives adjoint VEVs in the desired directions) 
can be constructed. And symmetries can be found that ensure its structure 
and the stabihty of the gauge hierarchy against possible Planck-scale effects. 
This has been shown in Ref. 16. 

It remains to show that symmetries can be found which render natural 
the full model, including both the Higgs and Yukawa parts of the superpo- 
tential. This is done in Appendix C. It should be emphasized that most of 
the technical difficulty of making the model natural has to do with the Higgs 
sector, and that the problem is not made significantly more difficult by the 
particular Yukawa structure assumed. 

6 Variants of the Root Model 

(a) A variant of the root model with e = — | 

As explained in Section 2(a), the root model consisting of the terms in 
M^pinor predicts that 

^ ^ Qi+Qi- 

There are two values of e which give the Georgi-Jarlskog result mP^/m'^ = 3, 
namely e = and e = — |. The first corresponds to Q = hn and gives the 
model proposed in Ref. 1 and studied in detail in previous sections. The 
second value gives an interesting variant that we shall now discuss briefly. 
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1 + 1^ 



(57) 



With e 
-(10) 



+ 6, \6\ -C 1, the mass matrices take the forms (see Eqs. (8) 



Uo = aTv 



Da = aTv' 







Usme 

5 






(-i + |5)sin^ 

{-I -IS) cose 





\ 

[I + IS) sine 
(I + §5) sine (1 + 15) cose J 

/ ' 



(58) 



(59) 



aTv' 



il- 15) sine 



(60) 






VO (-l + |5)sine (1 + IS) cose J 

There is one immediately apparent explanatory success: (m° / m°) / (m^ / m^) 
= 0(5). However, the contributions to Vcb from the up and down sectors no 



longer nearly cancel and ^ ^Jm^/m'^. The situation is thus the reverse 
of the model with |e| <^ 1. There the prediction for V^^ is good and that for 
{m^^/m^)/{mym^) is bad and has to be corrected by another mechanism. 
Here {m^^/mf) /{mym^) works well, but has to be corrected by some 
other mechanism. The question arises: How could V^^ be corrected without 
disturbing the Georgi-Jarlskog relation? One possibility is the following. 
All the mass matrices coming from PFspinor have the form 

/ 



mo 



g sin e 

V g^sin^ {Q^ + Q)cose 



(61) 



This means that 



while 



det23i^o Qi+Qi 



m° m^ml det23-Do Qd'^Qd' 



Qd" + Qd Qu" + Qu, 



i&ne. 



(62) 



(63) 
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Consider adding to each matrix a contribution 

/ 



AM = mo 



A ■ sin^ e A ■ sin ^ cos ^ | , (64) 
\0 A • sin 6* cos 6* A • cos^ 9 



with A the same at least for D and L. Then D33 and L33 remain equal, 
and so, for small ^, = m° remains vahd and ^ = still holds. 

' TOT dct23 D 

But dct23(Afo + AM) = (Asin2^)((Q^ + Q)cos^ + Acos^^) - (Qsin^ + 
A sin cos ^^) (Q'^ sin + A sin cos ^^) = — QQ'^ sin^ 9, which is unaffected by 
the addition of AM! Thus, the Georgi-Jarlskog relation remains good. On 
the other hand, for small 9, 

Qd<^ + A cos 9 Qu<^ + A cos 9 



tan^, (65) 



_{Qd- + Qd) + Acos9 {Que + Qu) + A cos 9 _ 

which is very greatly affected by the addition of AM, and can easily be made 
small. 

How can such a AM arise? In the original root model, defined by VFgpinor, 
the SU{2)i X f/(l)y-breaking term connects I]jajl6i = sin 6*162 + cos ^163 
with 16, which in turn mixes with I]j&il6j = I63. This gives the form of 
Mq. Clearly, if there were effectively a term of the form Y^ijO-iO-j^^i^^j it 
would give a contribution of the form AM. Such a term can arise in a mod- 
ified version of the root model. Imagine that instead of a simple vectorlike 
pair of family, there were two such pairs: 16 + 16 + 16' + 16'. Consider 
the superpotential (MT616 + Ei &iT616i45//) +{M'Tq'1& + Y.ia{iQ' lQ.il h) 
+glQ IQ'IQh + hlQ'lQ'lQu- Clearly the 16 and 16' mix with Ei ^il6j and 
Ei ail6j, respectively. Then the last two terms give, effectively, contributions 
of the form Mq and AM. 
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This variant root model has several features which make it seem less 
attractive than the model with e = 0. In the latter model the choice of small e 
explains three facts (the smallness of Vcb, the 2nd to 3rd generation hierarchy, 
and rnP^/rnPg = 3) and also helps explain the suppression of Higgsino- mediated 
proton decay. (In the long version it also plays a crucial role in explaining the 
suppression of m°/m°.) In the variant root model the choice of e explains only 
two facts (the smallness of m°/m° and m°/m° = 3). Secondly, the variant 
model is less economical by virtue of the introduction of the extra pair of 
family and anti-family, and involves a reintroduction of a family symmetry 
of sorts, since the 16' + 16' must be distinguished from the 16 + 16. Finally, 
the point e = — |, unlike the point e = 0, does not seem to correspond to a 
group-theoretically interesting direction. Nevertheless, perhaps some of the 
ideas in this section are capable of further interesting development. 

(b) A variant way to suppress m^/m^ 

In Sections 2(c) and (d) two ways to break the proportionality relation 
m°/m° = m°/m° based on a multiplicative correction to the mass matrices 
were described. They give the two versions of the model. Another possibility 
is a mechanism based on an additive correction to the mass matrices given 
in Eqs. (8) - (11). (We are now assuming again that e = 0.) Adding 
something to D to increase wPg/m^ without disturbing V^^ or the Georgi- 
Jarlskog relation seems difficult if not impossible. On the other hand, a 
small additive correction to U could approximately cancel off the small f/23 
element without significantly affecting U^2 or C/33. Then wP^/m^ would be 
suppressed without affecting V^^ or any of the relations that come from the 
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forms of D and L. 

There is only one difficulty with this idea: such an additive contribution 
to U would in general make U be rank 3 and, possibly, make too large. 
Since AC/23 = -[^0,23 = (m°m?/mg)/ sin ^ > m°, the form of AC/ must be 
such that Af/11 < AC/23. 

One way to ensure that AC/ has a form that maintains the rank 2 na- 
ture of Uq is by extending the root model (H^spinor) in a way analogous to 
what was described in Section 6(a). Introduce, as there, two pairs of fam- 
ily plus anti-family instead of one. Consider the superpotential W^^^^^j. — 
(AfT616 + &iT616i45^) +(M'T6'l6' + Ei Oiie'iea//) +^16 le'lOj^ + 
Ei /il6il6'16j^l6j^. The term ^flG 16'lOij gives effectively a contribution 
of the same form as Eq. (4). The last term contributes to U {ii {1{16h)) and 
(5(16ij)) are both non-zero), but not to D or L. The sum of the last two 
terms effectively involves only two linear combinations of 16's, and therefore 
the resulting total U is still rank 2. 

The first generation can still be given mass by adding IVvector- The present 
variant would dispense with need for VFadjoint; however. The main drawback 
compared to the long version of the model described in Section 2(d) is that the 
predictions of neutrino mixing angles are lost because of the extra parameters 

Appendix A 

Derivation of exact mass matrices 



38 



Consider the matrix 



where M and M' contain elements of order Mqut and m and mo contain 
only Weak scale entries. M. may be block-diagonalized as follows. 

^^^^^-[ (MtM + M'tM')^ J' ^^^^ 

where 

^ ( ^ 7 ^ (MtM + M'tM')-^ ) ( M'^ ) ' 

and 

/ -{mlM' + m^M){MW + M'^M'Y 

y [M^M + M'tM')-HM'tmo + M^m) I 

(69) 

Here x = M'M"^. Terms of order (Mweak/^^^GUT)^ have been dropped. This 
sec-saw form will now be applied to the case of the mass matrix of the charge- 
(— |) quarks in the short version of the model. 

We will define the superheavy linear combination of 16's to be 16' = 

. ^ , [16 + TQI63I, and the light orthogonal combination to be I60 = 

/ \ J -TQ16 + I63]. Then 

16 = \ [16'-rQ16l], 

I63 = - \ [+rQ16' + 16'3]. ^'^^ 

Define di, ^2, d'^, d', d'^, g, and g' to be the charge-(— |) quarks in the I61, 
I62, I63, 16', 16, 10, and 10', respectively; and define d^, ^2, d'^, d"^, d, g'^, 
and g"^ to be the charge-(|) antiquarks in the same representations. Then 
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by substituting Eq. (70) into Eqs. (1) and (2), and restricting attention to 
the down-type quarks, one gets the following 7x7 mass matrix 



{dl,d^„dld'^,d,g',g'') 



mo M' 
m M 



( di \ 

d2 

d' 
9 

\9' J 



where 



mo 



M' 













^ 











-aTv'QdSg/Nd 




]■ 




-aTv'QdcSe/Ndc 


-aTv'iQdc + Qd)ce/NdcNd j 






1 ° 


av'sg/Ndc 


av'{l-T^QacQd)ce/Nd.Nd \ 



















m — 


CiV' 


C2V' 






1 








c'^v'/Ndc 


/ 




f 







CiVr 


c'iVr 


\ 




av sg/Nd 


C2VR 


C2VR 



V av'{l - T^Qd'^Qd)ce/NdcNd c^VR/Ndc c'^vn/Ndc J 



(71) 



(72) 



(73) 



(74) 



/ aTv'iQdc + Qd)ce/NdcNd MNdc c^TQdcVn/Ndc c'^TQdcVn/Ndc \ 



M 



MNdc 
c,TQdv'/Nd 
c'^TQdi'/Nd 












+d{ibH) 





-d{ibH) 




V c'^TQdi'/Nd +d(45^) / 

(75) 

and where vr = (l(16i^)), v' = (5(16//)). v' = (5(10//)), cq = cos 6', and 
sg = sin^. Then by Eq. (57) 



Do = mo. 



(76) 
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and 

D ^ {I + xx^)-^mo- M'M-^m). (77) 
Comparing with Eq. (19) one sees that 

Arfc = xx^ = (M'M-^){M'M-y. (78) 

Multiplying out this expression using Eqs. (74) and (75) (one can neglect the 
terms of order v' and v' in M' and M) gives the result in Eqs. (22) and (23). 
And multiplying out the expression M'M~^m gives the flavor-antisymmetric 
piece in Eq. (13). 

The same kind of calculation gives the mixing matrices in the long version 
of the model as well. 

Appendix B 

Suppression of mc/rrit in long version of model 

Prom Eq. (27) one has 

A„c = E*E'^ + E'*E'^. (79) 

Define 

1 

2 

(80) 
So that 



p 


2 


E 


2 

+ 


E' 


2 

+ 


ExE' 



P = ExE' 



1 + 



E 



+ 



E' 



Ex E' 
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We will assume that 



E and E' are large compared to one (in fact, of order 
1/e), so that there is a hierarchy: P ~ E'"^ E, E' 1. Further, define 

(82) 



F = P* X E, 
F' = P* X E'. 



Then one can write an exact expression for the square of the mixing matrix 
(/ + A„c)-^: 

(7 + A„c)-i = YTiPf ^ ^ ^'^'^ ^ ^^^1 ■ ^^^^ 

This expression is easily checked by going to a particular basis. Since it is in 



a rotationally invariant form, it is true in any basis. 
One can write the full matrix U as 



with 



where 



U ^ {I + Auc)-Wo, 




Un=\ 7] I = (0,A,S) 
sin 9 cos 6 



rj = -esin^, 
5 



(84) 



(85) 



(86) 




Then the matrix Uq has eigenvalues 



mtfl 



-Tj sin 6 



AxB 



(87) 



(88) 



(89) 
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To find the eigenvalues of U one considers 



c/tc/ = uI{i + k^)-^Uq 

( \ 



1+IPI^ 



At 



(/ + FFt + F'F't + PPt) . (0, A, B) 



(90) 



Tfien 



1 



i+ipr 
1 

1 



It 


p 


2 

+ 


It 


P 


2 

+ 


St 


p 


2 

+ 


Pt 


■P 


2 

+ 



+ 



+ 



A 
B 



A^ -F' 
B^-F' 

^(At • PPt • P + At . ^^t . p + At • P'P't • P + At • P). 

(91) 

Tfie first row and column are zero in our approximation. One then has 

ml / ml q'^ ml ^ tT{WU) 



r-j 1 



l+IPJ 

~ i-Psr 



At • p + pt • p ) 



(92) 



i+|pr ' 

where we are keeping only the leading terms. And 



mlml = det23 WU. 



(93) 



Keeping only the leading terms (recalling the hierarchy P <^ F, F' <^ 1) one 
has after some algebra 

1 



mlml = 



At • Ppt • P - pt . PAt • P 



+ (P ^ F') 



(94) 



[i + \p\y 

The crucial point is that the leading terms in the brackets, which are O(P^), 
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cancel leaving O(P^). Further manipulations give 



(l+jPpj2 

AxB 



{A X By ■ (P X F) + (F F') 



{Ax By -E +{E ^ E') 



Or 



rricmt 



mc,omt,o 1 + |P| 



Eif + \E[f, 



and, since ^f- = Y^^p; suppression factor is given by 



mc,o/mt,o IP3 



Using the definition of P 

mjmt 



Ex E' 



mc,o/mt,o 



Eif + \E[\\ 



(95) 



(96) 



(97) 



(98) 



{E X E% 

This is the desired result. Notice that for large E and E' this goes as 
1/ min(£', F'), unless E x E' happens to point nearly along the 3 direction. 
In other words, generically the ratio rric/mt is suppressed. Since E, and E' 
are of order ^, rric/mt is suppressed by a factor of 0(e). 

Cases of A = rank 1 or rank 3 

The curious fact that mc/rrit is generally suppressed is a result of the fact 
that is rank 2. We can show this by doing the analogous calculation for 
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the cases where A is rank 1 and rank 3. If A is rank 1 it can be written 



A = E*E'^. 



(99) 



Then 



(/H-A)-i = / 



* 



l + \E\ 



,E*E 



(100) 



ParalleUing the calculation for the rank-2 case closely one finds 





El 


2^ 


/l + 


E 


mc,omt,o 1 + 


E, 




E2 


2 



(101) 



This is not small unless E happens to be nearly in the 2 direction, which 
agrees with intuitive expectation. 

If A is rank 3, then (7 + A)"^ is just an arbitrary 3-by-3 matrix which 
we will call M. A calculation similar to the preceding two gives 



rrir 



mc,omt,o 

Again, this is not small except for particular special choices of A. 



(102) 



A calculation of V^f^ 

Let the transformation of the left-handed charge-(|) quarks required to 
diagonalize U he 



/ 1 \ 
Vu^ c s* , (103) 
\ -s c J 

and the matrix required to diagonalize Uq be the same with s ^ sq and 
c — > cq. Similarly, let the transformation of the left-handed charge-^ quarks 
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required to diagonalize D be of the same form with s ^ s' and c — >■ c'. We 
take c, c', and cq to be real and s, s', and sq to be complex. Then 



- sc' - s'c 



(104) 



The primed quantities arc easy to find: s'/c' = D32/D33 = (1 — |e) tan^^ as 
can be seen from Eq. (15). (We ignore the effects of C23 and the Nf.) This 
gives 

1 _ 1 

(105) 



s' ^ (sin^(l- -eRCOs'^e),--eisine). 
5 5 



And So = {Uq)^2/ (t^)33 which is given by the same expressions with e — > — e, 
as can be seen from Eq. (14). To find s we can use the expressions derived 
for U'^U earlier in this Appendix: 



2|s|c 
— \s\ 



32 



(WU)U-{WU) 



(106) 
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args = arg(^7^^7)32. (107) 

This gives, after multiplying out the expressions for the elements of WU 
given above, 

O I ol ^ O I I I Ti^T^^ 



and 



where 



2 


s 


c 


2 




1 (co + ReK) 


c2- 




s 


2 ~ 2 
co- 




So 


^ + 2coReis: 



arg s — arg sq — ImK/ cq. 



P3P2* + -P^3i^2 + -F^^r 

\P3f + iFsf + \F'\^ + 1 



(108) 



(109) 



(110) 
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(Ill) 



These expressions can be combined to give after straightforward algebra 



Appendix C 

In order to find a symmetry that makes the superpotential of Eq. (56) natural 
and stable it is necessary to introduce some singlet superfields. This can be 
shown as follows. With the terms tr(Ai)^ and (>li)^, Ai cannot transform 
non-trivially except under a Z2. (Similarly for A2 and A^.) But it is crucial 
for the doublet-triplet splitting that (T2)^ be allowed, while (Ti)^ is forbidden, 
which implies a non-trivial relative transformation of these two fields and 
hence of Ai because of the presence of the term T1A1T2. Thus the terms 
tr(Ai)^ and (^1)^ must be replaced with some form that allows a non-trivial 
transformation of Ai. The simplest possibility is to insert singlet superfields 
with non-trivial transformation properties. For example, 0itr(Ai)'^/Mpj. For 
similar reasons it is convenient to introduce singlets into other terms as well. 

There is no point in trying to find the most elegant or simplest combina- 
tion of symmetries that works. Rather here it will only be shown that some 




sin ^ cos ^(-e — K/ cos9), 



(112) 



or 




(113) 
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symmetry can be found. The easiest way to do this is to restrict the search 
to a single U{1) symmetry, and to introduce singlet fields where convenient 
to make a needed term allowed. No attempt has been made to economize on 
these singlets. 

Consider, then, the following set of fields. Its U{1) charge is given in 
parentheses after the name of the field. Adjoint Higgs: Ai{—a2 — a3), ^2(^2), 
^3(03); Fundamental Higgs: Ti(ii), T2{t2; Spinor Higgs: C(— c), C{c — x); 
Singlet Higgs: P{-ti - ^2 + 02 + 03), Q(-2t2/3), i?(r), S{s), X{x), 
Matter Spinors: Fi{e), F{—e — ti), F{e -\- — s); Matter Fundamentals: 
T{-a^/2 - y/2), T'{-a^/2 - y/2); Matter Adjoints: ^(-02/2), ^'(-02/2). 
These charges are not all independent, but satisfy c = (03 — 02 + 2a; + |/)/4, 
and ti — —03/2 — 802/2 + X — y/2 + s. With these charge assignments the 
following terms are allowed: 

W = tr(A)Vi/M2i+(A)Vi 

+ tr(A2)V2/M2j+ (^2)202 

+ XCC + XCC(^303)V^P1 + (^303)7^3 

+ iiiAiA^A^) , . 

+ Ti4iT2(P/Mpi) + (T2)2(gVM|i) + riCC(i?/Mpi) ^ ' 
+ [s{FF) + {FF,)A2 + (FF,)Ti 

+F(Tr')-l3/Mpi t_(TF,)C + (T'F,)C 

+{AA')A2 + {AFi)C + (A'F,)C]. 

All the terms have coefficients of order unity that are not shown. The VEVs 
of the fields are of the following orders of magnitude. (^1), (^2), (C), (C) ~ 
Mgut, W ~ Mpi, (X),(y),(P),(g),(i?),(5) ~ Mgut, {(l>i) ~ Mpi, and 
~ Mgut/Mpi. (Thus I/M3 must be of order M^ut/-^pi-) 
The most dangerous terms for the hierarchy are those that lead effec- 
tively to TiCC or to linear terms for Ai. The lowest term involving TiCC 
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is TiCC[A203-RF]/Mpi (assuming a field Y with opposite quantum numbers 
to Y exists). This gives effectively {Mqut / Mpi)^TiC C , which is sufficiently 
suppressed. The lowest dimension terms linear in Ai (that would desta- 
bilize (Ai); note that ti{AiA2A3) does not) are Ai{XCC){RSY) /M^^ and 
{AiA2As){XCC)/Mpi. The first is harmless, but the second gives effectively 
(Mqut/^pi)^i- Since the mass of Ai is necessarily of order M^^jr^/Mpi, one 
has 6{Ai)/{Ai) ~ Mqut/^pi- What is needed is 10"'', so that this danger- 
ous term must be assumed to have a dimensionless coefficient that is of order 
10~^. All other potentially dangerous terms are sufficiently suppressed if 
one assumes that all terms are suppressed only by dimensionally appropriate 
powers of the Planck mass. 
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Table I: 



Matter : 




F 


F 


T 


r 


A 


A 


Higgs : 


A2 


^3 


Ti 


C 


C 


^3 


1 




z 


z' 


z' 


z 


z 




z 




z' 


z 


z' 



Figure Captions 

Fig.l: The Wgpinor contribution to light fermion masses. 
Fig. 2: The Wvector contribution to hght fermion masses. 
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